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We study the liquid-solid transition in a collection of inter- 
acting particles moving through a dissipative medium under 
the action of a constant, spatially uniform external force, e.g. 
a charge-stabilized suspension in a fluidized bed or a flux- 
point lattice moving through a thin, current-carrying slab of 
type II superconductor. The mobility of a given region in 
these systems is in general a function of the local concentra- 
tion. We show that the structure factor peak is suppressed in 
an anisotropic manner as a result of this effect, resulting in a 
shift of the crystal-liquid phase boundary towards the crystal 
side. A nonequilibrium phase diagram is presented. 

05.40+j, 64.70.Dv, 82.70.Dd, 47.15.Gf, 



I. INTRODUCTION AND RESULTS 

How a liquid freezes to form a solid is a problem of 
longstanding interest to physicists. Freezing at thermal 
equilibrium can now be said to be quite well understood 
in terms of a first-principles order-parameter theory [Q. 
The effect of steady nonequilibrium driving on the freez- 
ing transition has been studied in the context of sheared 
colloids A uniformly driven state, with no gradi- 

ents on the average, as would arise in steady sedimenta- 
tion j6|,|7j or the fluidized bed |]§|-|l0| is conceptually sim- 
pler: unlike in shear flow, an order-parameter description 
of the mean state of the driven crystal is straightforward. 
Nonetheless, we know of no studies, experimental or the- 
oretical, of the effect of sedimentation on the liquid-solid 
phase boundary in a suspension settling uniformly under 
gravity, or of the freezing of moving flux lattices in the 
absence of quenched disorder. 

A study of freezing in the presence of uniform driv- 
ing is also of interest as a problem in the general class 
of driven diffusive systems. Most of the investigations 
in this field [llj| are either away from any equilibrium 
phase transition or in the neighbourhood of a second or- 
der phase transition. The effect of uniform driving on 
a first order phase transition has received very little at- 
tention [p^[ ; the freezing transition, in particular holds 
interesting possibilities since the ordering is at non-zero 
wavenumber. 

In this paper we consider a collection of interact- 
ing particles in a liquid state, i.e., with strong short- 
range translational correlations, and ask what happens 



to the liquid-solid transition in the presence of an im- 
posed steady mean translational motion of the particles 
with respect to the ambient medium. Our model should 
apply not only to charge-stabilized colloidal suspensions 
p3| , but also to the effect of a steady current on the 
freezing transition of flux-point liquids Jl4[ in thin slabs 
of clean type II superconductors. The main ingredient in 
our theory is the concentration dependence of the mobil- 
ity of particles in a suspension, or of the flux-points in 
a superconductor. In colloidal suspensions, this depen- 
dence arises as a result of the hydrodynamic interactions 
between particles, while in the flux-liquid, it is the elec- 
tromagnetic interaction that is primarily responsible for 
the effect. Although our approach is fairly general, we 
use parameter values appropriate to the specific case of 
sedimenting colloidal suspensions in the rest of the paper. 

Our main results are as follows: As the mean concen- 
tration of the suspension is increased, the hydrodynamic 
interaction causes it to settle more slowly p, ^5|Jlr| | . We 
show that this leads to a nonlincarity in the equations of 
motion whose effect is to reduce the height of the struc- 
ture factor relative to its equilibrium value. The suppres- 
sion is strongly anisotropic, largest for wavevectors along 
the direction of mean drift, and vanishing for wavevec- 
tors normal to it. The reduction is most pronounced at 
the location of the peak, which means the strong short- 
range correlations in the liquid, which are the precursor 
of translational order, are suppressed. This inevitably 
leads to a shift of the phase boundary towards the crystal 
side (i.e., favours the liquid over the crystalline phase), 
by an amount proportional to the square of the Peclet 
number, in a manner reminiscent of a proposed mecha- 
nism Q for shear-induced melting. This nonequilibrium 
phase transition now assumes additional importance in 
the wake of recent experiments on charged fluidized beds 
|p[- ^0| , ^5[ where it should be possible to test our predic- 
tions. Constraints on the types of systems to which our 
theor y should apply are discussed at the end of section 

vn| 

The paper is organised as follows. Section || gives an 
overview of equilibrium freezing. The equations describ- 
ing the motion of a sedimenting suspension are derived 
in section [II. In section IV we discuss in a general set- 
ting the prescription used in this paper to describe a 
nonequilibrium phase transition. The perturbative cal- 
culation leading to the equal time correlation function 
in the driven system is the content of section |v|, and 
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our numerical results f or th e structure factor are given in 
section VI. In section VII we derive our results for the 



freezing transition, and conclude in section VIII 



II. EQUILIBRIUM FREEZING 

The equilibrium phase-transition from a liquid to a 
crystal is the spontaneous appearance of a static, spa- 
tially periodic modulation in the density field of the liq- 
uid. An order parameter theory of this transition p"7| , |l8[ . 
is obtained by writing down the free-energy cost of such 
a modulation as a functional of the density field. 

The free energy cost of an inhomogeneous density field 
p(x) relative to the uniform liquid state of density po is 
given by 



F = F/k B T 

= \J dx C (|x-x'|)(p(x) 
P(x) 



rfx 



p(x)ln 



P(0) 



(p(x) - po) 



Po) 



(2.1) 



where c(r) is the direct co rrel ation function of the liq- 
uid. The second term in (2.1) is an ideal gas entropy, 
whereas the first comes from correlations in the interact- 
ing system. Expanding the logarithm and stopping at 
the second order in the density fluctuation <f> = p — po, 
gives the quadratic part of the free-energy in terms of the 
Fourier modes <j>(k) as: 



F 



dk 



1 



-I0(k) 



{2n f /OoSb(fc) 1 
in terms of the dimensionless structure factor 

So (*0 ' 



1 - p c(k) 



(2.2) 



(2.3) 



The density correlation function at equilibrium is related 
to the structure factor by the relation 



(0(k)0(-k)> = VpoSo(k) 



(2.4) 



where V is the system volume. 

All information about temperature and mean density 
is encoded in the liquid-state direct correlation function 
c(r). For a given c(r), p{r) corresponding to the thermo- 
dynamically stable state, within mean-field theory, is the 
global minimum of F[p], obtained in principle by solv- 
ing ^ = Q, and choosing that solution for which F is 
the smallest. At low densities and high temperatures, 
the global minimum is p(r) = po. For sufficiently dense, 
cold systems, this uniform liquid state is unstable rela- 
tive to a spatially periodic p(r), the crystal. Terms in F 
of order higher than quadratic in cf>(r) = p(r) — po are 
clearly essential to give such a solution. Although the 
best theories of freezing go the whole log, we shall 



sacrifice accuracy in favour of simplicity, and use a quar- 
tic polynomial approximation to F. Since we are looking 
for periodic density waves, it is convenient to work with 
{</>g}, the weights of the Bragg peaks at the reciprocal 
lattice vectors G, defined by 



0(k) = ^ G <5(k-G). 



(2.5) 



Since the peak of the liquid-state structure factor at its 
first maximum G — 2n/a is much higher than the subse- 
quent peaks, we restrict ourselves to Fourier amplitudes 
of 4> close to this peak. This is equivalent to consider- 
ing only the first shell of reciprocal lattice vectors for the 
crystal. Only the form of the S(k) near the peak matters 
then, and S~ 1 (G) can be replaced by a Lorentzian 



S-\G) = a + a(G 2 - fc 2 ) 2 
so that the peak height is 

S(k )=a-\ 



(2.6) 



(2.7) 



The Landau Free energy per particle is then given by 
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G1...G4 



Po 



3Gi+G 2 +G 3 ,0 

<5gi+g 2 +g 3 +g 4 ,o- 
(2.8) 



If aGo ^> a, then free energy costs strongly discour- 
age the Fourier amplitudes to stray from the peak of the 
structure factor, i.e |G| = fco- The free energy can then 
be expressed as a function of a single scalar order param- 
eter 



r = n 1/2 (pa/po) 



(2.9) 



where n is the number of reciprocal lattice vectors in the 
set (i.e. 6 for a triangular lattice, 12 for bcc etc). Note 
that at nonzero wavevector G, 0g and pg are identical. 
The free energy per particle per ksT is thus a polynomial 



F 



Nk B T 



1 2 
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with 



b = 2pn 



-1/2 



Po 2 B 



c = 6 



1 + 1 



Po 3 C. 



(2.10) 



(2.11a) 
(2.11b) 



Here p and q are numbers that depend on the type of lat- 
tice being considered: p is the number of such triplets or 
triangles that each RLV belongs to, and q is the number 
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of non-planar diamonds that the RLVs can form among 
themselves fis|| . The value of a at which the transition 
occurs is found by demanding that the free energy be 
zero at the crystalline minimum. This leads to the value 
of a at the transition, 



Pe = 



3c' 



(2.12) 



The cubic term favours lattice structures in which the 
RLVs form triangles in sets of three. The candidate lat- 
tices, within the one order parameter theory, are the ones 
whose RLVs correspond to edges of regular polyhedra 
or polygons with triangular faces, i.e. planar triangular 
(whose RLVs form a pair of triangles), face centred cubic 
(octahedron) and icosahedral (icosahedron) lattices. The 
theory predicts that all crystals should be triangular lat- 
tices in two dimensions and body centred cubic (which 
has a face centred cubic reciprocal lattice) in three di- 
mensions. Given the lack of sophistication, the theory is 
reasonably successful as almost all two dimensional crys- 
tals are triangular lattices, and many, if not all elemental 
solids in three dimensions are indeed BCC, at least close 
to freezing Jl7[ . 

The theory discussed above has several drawbacks: it 
neglects all structure in S(k) except that at the peak 
and assumes that all wavevectors have the same order 
parameter amplitude. These assumptions can be re- 
laxed and a more complete theory written down, but 
the Landau expansion itself is questionable, in fact in- 
accurate, since the order parameter jump at freezing is 
large. It would doubtless be better to generalize the first- 
principles density-wave approach of Ramakrishnan and 
Yussouff Q] to include the effects of sedimentation. In 
the present work, however, we have attempted only a 
nonequilibrium generalizations of the simple and intu- 
itive Landau-Alexander-McTague theory, as a first step 
in taking into account the effect of sedimentation on 
freezing. An extension of our approach to more realis- 
tic freezing theories should be fairly straightforward. 



III. EQUATIONS OF MOTION 

The frictional force on a spherical particle of radius R 
moving slowly in a fluid with shear viscosity rj with a 
uniform speed v is given by J2(| 



(3.1) 



The dimensionless parameter which gives a measure of 
the relative importance of sedimentation with respect to 
Brownian motion is the Peclet number, defined as the 
ratio of the time taken by the particle to diffuse a char- 
acteristic distance t to the time taken to settle the same 
distance under gravity, without diffusing. In terms of 
the mean settling speed ^stokes of the particle, the Peclet 
number is given by 



^Stokes I 
D ' 



(3.2) 



The appropriate dimensionless measure (which we will 
continue to call the Peclet number) of the relative im- 
portance of drift versus diffusion for this problem, will 
turn out to be proportional to the Peclet number defined 
at the scale of an interparticle spacing, but will incorpo- 
rate certain collective effects whose nature will become 
clear below after we construct our equations of motion. 

Although the single particle problem was solved in 
1851 ^0|, the sedimentation speed of a suspension was 
theoretically obtained as recently as 1972 |2l|j , and that 
too in the dilute limit. The speed of a strictly random 
suspension of volume fraction <fi turns out to be less than 
that of a single sphere: 



«*=(!- 6.550 + 0(4> 2 ))v st , 



okes 



(3.3) 



Thus, at a volume fraction of 10%, the settling speed is 
reduced to 35% of the Stokes value. It has been also 
observed experimentally jl5|,[l6) that the mean settling 
speed of a suspension is less than that of a single par- 
ticle, and a decreasing function of the concentration of 
suspended particles. In other words, the mobility of a 
region is a function of the local concentration. 

The sedimenting suspension that we speak of is an ide- 
alized bottomless tube of infinite height containing a sol- 
vent in which there is a steady downward flow of particles 
(denser than the solvent). An 'infinite' tube is realized 
in real life by means of a setup known as a fluidized bed 
|§- |l0| , ^5| , a vertical tube through which fluid is pushed 
up at a steady rate by an externally maintained pressure 
head, with particles suspended stably on the average by 
the balance between gravity and the viscous drag of the 
upfiow of fluid. The problem of the steadily sedimenting 
particle in the comoving frame of reference is identical, 
in the infinite container limit, or far from the walls of the 
container, to the fluidized bed problem in the lab frame. 
For very dilute solutions, the reduced weight of the parti- 
cle (i.e. the weight minus the buoyant force) is balanced 
by the single particle Stokes friction. At higher concen- 
trations, hydrodynamic interactions between the spheres 
p2| , p3"t becomes important, leading to a hindered settling 
speed, the effect of which we investigate. 

Consider now a collection of interacting particles, on 
the average in a liquid state with uniform number density 
Po, moving through a frictional medium. In the absence 
of external forces (such as gravity) the number density 
p(x) follows the equation: 



dp 2 SF 



v, 



(3.4) 



where Tq is the collective diffusion constant, r\ is a Gaus- 
sian white noise consistent with thermal equilibrium: 

(?7(x, t)?7(x', t')) = 2fc B Tr V 2 ( 5(x - x.')S{t - t'), (3.5) 



and F is the free-energy functional defined in 
equation is simply the equation of continuity 



This 



3 



dp 

at 



= -v.j 



(3.6) 



with a current with two pieces: a mean (J) = —ToVip 
and a fluctuating piece f , uncorrelated in space and time. 
The chemical potential if) is related to the Gibbs free en- 
ergy F by the thermodynamic relation tp = ^£ which 



gives equation (3.4) with a conse rvin g noise r; = — V.f. 
Note that the e quat ion of motion 3.4 combined with the 
noise statistics ( |3.5[ ) will take the system to an equilib- 
rium state with a probability distribution 



P 



-F/k B T 



(3.7) 



In the presence of a constant external force, which we 
take for specificity to be a uniform gravitational field 
g = gz, the particles sediment with a steady mean drift 
velocity v = /i(p)mg |24|) , where p is the concentration- 
dependent mobility, which in the low density limit, is 
expected to go to the Stokes value 6 ^ vR - The current 
then acquires an extra piece pv leading to an additional 
term of the form 



(3.8) 



on the right hand side of equation (3.4). Here mg is the 
reduced or Archimedean weight of the solute particles, 



mg = m g - m displaced g, 



(3.9) 



with m being the mass of the particle, and ^displaced 
the mass of solvent displaced by submerging it. 

Let us now consider concentration fluctuations in the 
suspension around the mean density po : P = Pa + 6 p. 
Taylor expanding the mobility about the mean value po 
to second order, we find the equation of motion for the 
density field to be 



0_ 

dt 



Sp + mg(p Q + po 



dp 
dp 



d 



Sp 



= r v - — P2i^-Sp 

op oz 



(3.10) 



In the rest of this paper, we will do all our calculati ons in 
this reference frame, so that the left hand side of ( 3.10| ) 
can be replaced by a time derivative. Note that this frame 
is not the frame in which the mean velocity of the parti- 
cles is zero. For a colloidal suspension, since p decreases 
on increasing concentration, the frame moves with a ve- 
locity less than the mean velocity. Thus one would see 
a net drift along gravity in this frame. However, the 
quantities of interest to us, i.e., equal time correlation 
functions, are invariant under a Galilean transformation, 
so it doesn't matter in which frame we do our calculation 

Before beginning the calculations, we change variables 
to the non-dimensionalized fluctuation of the concentra- 
tion from the mean, 6 = p ~ Po , and the dimcnsionless 

P° 

free energy F = F/ksT, in terms of which the equation 
of motion in the comoving frame becomes 



d(f> KraSF , d , a 

i = DV H +x d-z* +r >> 

(?7(x, f)r/(x', t)) = -L>V 2 (5(x - x)5(t - t'), 



where 



D 



Po 



(3.13a) 
(3.13b) 

(3.14) 



is the collective diffusion constant, and 



-Pomg 



1 ( E 



d I In 



(3.15) 



where 



is a velocity scale set by the Stokes velocity modified by 
the concentration-dependence of the mobility. Note that 
the piece pomg outside the brackets is the mean settling 
speed. The term within the brackets is a dimensionless 
number that vanishes in the dilute limit. For sediment- 
ing suspensions, the mobility is a decreasing function of 
concentration, so A > 0. Since D and A are a diffusion 
constant and a velocity respectively, together they give 
the appropriate dimensionless coupling (which we will 
call the Peclet number Pe) for the problem: 



P2 = mg 



Pomg 
Po 



dp 
~d~ P 



Po 



a i 



i>t) 



d 2 p 

dp 2 



(3.11) 



The left hand side of (3.10) is simply the Eulerian deriva- 



tive for a reference frame moving with a velocity 



v fr 



g(po + ( 



d(lnp) , 
dp 



d 



(Sp) z. 



(3.12) 



D 



ppmglptf 
ksTTo 







(3.16a) 
(3.16b) 



where £ is the average interparticle spacing in the sus- 
pension. Note that we are using I rather than the particle 
size as the length scale to define the Peclet number. This 
is reasonable since we are considering collective rather 
than single particle effects. In terms of the single par- 
ticle diffusion constant Dq, which is related to the bare 
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mobility by the Einstein relation Dq 
number can be written as 



ksTfio, the Peclet 



p e =ra 

D 



d I 



( mgR \ 
\k B T 



(3.17) 



The last piece (in parentheses) in ( pM7| ) is the particle- 
size Peclet number for a single sphere, which is simply 
the ratio of the gravitational potential energy drop across 
the size of the particle to the thermal energy k^T . The 
other three terms are dimensionless numbers that bring 
in the collective effects. 

The hydrod ynam ic interactions appear through the 
third piece in 3.17 (in square brackets). It should be 
noted that this term is indeed nonzero even if the mobil- 
ity is a linear function of concentration. 

We pause here to do some power-counting. The 
change in mobility due to a number density fluctua- 
tion depends on the change in volume fraction. The 
Taylor expansion of mobility, for example, is given by 
p = fio [l — ctR d (p — po) + . . .1 in d dimensions, where 
a is a dimensionless number that can be found from ex- 
periments. To the lowest order in the polyb all size R, 
therefore, the piece in square brackets in 3.17 is propor- 
tional to R d . The Peclet number thus goes as mR d for 
particle of mass m and radius R. As R — ► 0, the hydro- 
dynamic interaction between the spheres goes away, and 
the Peclet number goes to zero as R d if we keep the mass 
fixed, and as R 2d if we change the radius of the particle 
keeping the density fixed. 



The equation of motion ( 3.13 ) of the density fluctua- 
tion is the equilibrium dynamical equation with an ad- 
ditional nonlinearity of the Burgers type p6| . Qualita- 
tively, the nonlinear term may be understood in terms of 
a lattice gas model [Q , where particles are driven in one 
direction by a constant field, but a particle is allowed to 
hop to the next site only if it is empty (no double occu- 
pancy). This microscopically models a mobility depend- 
ing on the local density - in regions where the density is 
large, the probability of hopping being obstructed by the 
presence of other particles is larger, and hence the mo- 
bility smaller. The current at site n is then proportional 
to p n (l — p n +i), where p n is the discretized density in a 
bin at n. On coarse graining this will give two terms: a 
term ~ <f>, which can be removed by a Galilean transfor- 
mation; and another ~ — 2 , which persists even in the 
comoving frame. 



Before moving on to use (3.13) to calculate the corre- 
lations in the suspension, let us briefly discuss the limi- 
tations of our model. First, we have only considered one 
effect of sedimentation, namely the dependence of con- 
centration on mobility. This is a coarse-grained descrip- 
tion of the hydrodynamic interactions, reliable only at 
long wavelengths. A better theory should consider, for 
example, gradient-dependence in the mobility, i.e., the 
diffusivity at wavenumber k should have corrections of 
O(kR) for a suspension of particles of size R. Since we 



are interested in ordering phenomena at finite wavevec- 
tors ~ fcoj our theory should give reliable results for 
koR << 1 which is satisfied, for example, for strongly 
charge stabilized colloids where the ordering length scale 
fco 1 is of order the screening length, £ 3> R. In addition, 
we have ignored the dependence of the mobility on more 
detailed aspects of the local distribution of particles, e.g., 
its anisotropy, a point which we comment on at the end 
of this paper. 

Secondly, we have also ignored the concentration- 
dependence of the diffusion constant D in (3.13a). 
Strictly speaking, D should have the same concentra- 
tion dependence as po, since the two are related by the 
Einstein formula. A correct formulation of the prob- 
lem with D{4>) would require a multiplicative noise (of 



the form V.[y D(4>) f] where f is a spatiotemporally 

white noise with unit variance) to make the equations 
of motion consistent with equilibrium. This problem is 
technically much more difficult [^,^8), and we believe 
it shouldn't change the results dramatically for the fol- 
lowing reason. A concentration-dependent diffusion con- 
stant coupled with a consistent noise of the above type 
is an equilibrium effect, and an (effective) concentration 
independent D with ordinary Gaussian noise leads to 
the same equilibrium distribution. The nonlinear term 
in (3.13a) is the only driving term arising from concen- 
tration dependence. In any case, ( J3. 13 ) is the simplest 
driven-diffusive model with a liquid-solid transition and 
the correct equilibrium limit, and thus merits study in 
its own right. 

Lastly, we have also ignored the long-ranged effects 
of the hydrodynamic interaction. Our description corre- 
sponds to a quasi two-dimensional suspension confined 
between parallel walls with spacing h. The walls cut off 
the hydrodynamic forces on lateral scales much larger 
than h, making them effectively short range. The local 
effects of the hydrodynamic interactions are, however, 
accounted for by the concentration dependent mobility. 
To describe the three dimensional problem we have to 
explicitly include the the hydrodynamic velocity-field of 
the suspension in the description. Although this could 
have far-reaching effects at small wavenumber, it should 
be relatively unimportant for the statics of crystalliza- 
tion, whose physics is dominated by wavenumbers near 
the structure-factor peak. 



IV. PHASE TRANSITIONS IN A 
NONEQUILIBRIUM PROBLEM 

Since we are interested in the freezing transition, we 
set up the equation of motion (3.13) for the density 
field, with a free energy functional that leads to a freez- 
ing transition in equilibrium. The simplest of these is 
the Landau-Alexander-McTague free energy functional 
which we have discussed in section ||. A spatially inho- 
mogeneous density fluctuation <fi(r) (Fourier transform 
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k)) has, relative to the uniform fluid, a free energy 
dk 

dp 

"J {2n) d J (I 

1 



F[0(k)] = i / ^ 1 ^(k)oi-k) 
2J (2tt) PoAo(fc) 
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(27T) 



rfp 



-k)0(p)0(k - p) 
dp' 



(27r) a 7 (27r) a 7 (27r) a 
(-k)0(p)0(p')0(k - p - p') (4.1) 



where So(k) is the equilibrium structure factor of the sus- 
pension. The density correlation function at equilibrium 
calculated using the above free energy is, to quadratic 
order, 



(0(ki)0(k a )> = Po S (k)S(k 1 



(4.2) 



The equation of motion ( |3.13| ) becomes in reciprocal 
space 



dp 



(4.3a) 



+i—k z I -^0(p)0(k-p)+?/(k), 
(r/(k, t)r?(k', *')) = (2ir) d 2Dk 2 p 6(k + k')S(t - t 1 ). (4.3b) 



Using t he fo rm ( |4.lD for the free energy, the equation of 
motion (4.3a) becomes 



f) f)k 2 A f dn 

-m = - Tm m + r-k z / p <^(p)0( k - P 

ot 6o(KJ po 7 (2tt) 



-,7(k). 



(4.4) 



where the contribution of anharmonic terms in the free 
energy have not been put in. This is because we are 
itcrcstcd in the corrections to the equilibrium bare cor- 
relation funtions due to the driving term. 

It is easy to sec that the additional nonlinear term in 
equation ( 3.1 3| ) or (4.4) cannot be expressed in a form 
so as to be absorbed into an additional term in the free 
energy. This is done by checking that the functional curl 
of the term is nonvanishing. We will call this term in 
(4.4) the sedimentation nonlincarity. 

To describe phase transitions in this nonequilibrium 
problem, we follow earlier work on critical phenomena in 
the presence of shear [^9| and shear-induced melting of 
colloidal crystals 0. We assume that the Fokker-Planck 
equa tion that derives from the above Langevin equation 
(4.4) has a steady solution in the long time limit, which 
we write as 



P[4>] ~ e" rM (4.5) 
and expand the functional T in a Taylor expansion 



T[cf>} = ~r 2 (l, 2)0(1)0(2) + Ir 3 (l, 2, 3)0(1)0(2)0(3) + . . 

(4.6) 

where 1,2 etc. are arguments, for example spatial po- 
sitions Xi,X2 etc. The summation convention is implied 
for repeated arguments. The T n are the vertex functions. 
The calculation of moments from this probability func- 
tional gives, for the connected parts, the relations |3(| 

(0(1)0(2)) = r 2 (l,2r\ (4.7a) 

(0(1)0(2)0(3)) =-r 3 (l',2',3')(0(l)0(l')> x 

(0(2)0(2')) (0(3)0(3')) , (4.7b) 

and so on, which are the usual relations between the 
connected correlation functions and vertex parts. The 
only difference here is that the steady state probability 
distribution is not the inverse exponential of the equi- 
librium free energy functional. The vertex functions I\ 
above have to be therefore calculated from the steady 
state probability distribution P. However, the correla- 
tion functions can be calculated from the Langevin equa- 
tion itself by demanding that the time deri vativ e of all 
moments is zero, and these may be used in (17) to cal- 
culate the vertex functions upto any desired order, which 
gives the steady state P to the same order in the Taylor 
expansion. For instance, 

-l 



Wki) 
Wki) 



(k 2 )) = r 2 (ki,k a ) 
(k 2 )0(k 3 )) = 



(4.8a) 



f dki f dk 2 ' f dk 3 ' , , , 

(0(k 1 ')0(k 1 )) (0(k 2 ')0(k 2 )) (0(k 3 ')0(k 3 )) . (4.8b) 

The second of these relations can be rewritten in terms 
of the structure factor as 

(0(kx)0(k 2 )0(k 3 )) = -r 3 (k 1 ,k 2 ,k 3 ) Po 3 5(fc 1 )5(fc 2 )5(fc 3 ). 

(4.9) 

At thermal equilibrium, i.e., in the absence of driving, 
the vertex functions calculated to tree level are simply 
the inverse correlation function and the higher vertices 
in the free energy: 

' ^(ki + ka), 



r 2 (ki,k a ) 

r 3 (ki,k 2 ,k 3 ) 
r 4 (ki,k 2 ,k 3 ,k 4 ) 



poSo(ki) 
= -b5(ki +k 2 +k 3 ), 
= c5(ki +k 2 + k 3 +k 4 ). 



(4.10a) 

(4.10b) 
(4.10c) 



We calculate the perturbative correction to the vertices 
due to the sedimentation nonlincarity, from w hich we ob- 
tain the steady state probability distribution (4.5). This 
P can be treated as a mean-field probability distribution 
which may be maximized to find the stable phase at any 
given temperature. In other words, T[0] can be treated 
as a nonequilibrium free energy. Thus the perturbative 
corrections to the vertices can be interpreted as addi- 
tional terms generated by sedimentation in an effective 
free energy. 
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V. PERTURBATION THEORY 



We start with the linearized equation of motion (4.4) 
for the density field with only the sedimentation nonlin- 
earity, 



0_ 

at 



Dk 2 

. A , 

Po 



0(k) 



dp 



(p)0(k-p)+ry(k), (5.1) 



and calculate perturbative correction to the v ertex func- 
tions. The noise statistics is as given in ( 4.3b| ) which, for 
momentum labels k and — k reads 



(ry(k, t)r)(-k, t')) = 2Dk 2 Po VS(t - t') 



(5.2) 



where V is the volume of the system. To the lowest (i.e. 
first) order, there is a correction to the three-point ver- 
tex, whereas the structure factor gets modified at second 
order in Peclet number. The four-point vertex T4 will 
also have a correction at second order, but we ignore this 
correction as this is simply a stabilizing term and small 
changes in it will not affect the physics of freezing signif- 
icantly. 

In this section, we set up a hierarchy of equations for 
the moments of tf>, which we close by decoupling the four- 
point function into a product of two two-point functions 
( Hartree approximation). Multiplying the equation of 
motion (4.4) by 0(k), symmetrizing in fc, —k and setting 
time derivative of the two-point function to zero, we ob- 
tain for equal time a relation between the two-point and 
the three-point functions: 

f)k 2 

2^r- <#-k)0(k)) = 



S (k) 



l ~ k * I 77^K0(p)0(k - P )<K-k)> - (k - -k)] 
Po J (2tt) 

(5.3) 



(2nT 
2NDk 2 , 



where we have used the equal-time correlation between 
the noise and the density fluctuation, 



(<P(-k,t)r ] (-k > t))=NDk 2 , 



(5.4) 



which is proved in detail in the appendix. Here N = poV 
is the total number of particles in the system. Since the 
correlatio n fu nction is symmetric under interchange of k 
and k, (5.3) can be rewritten as 



Dk 2 
S (k) 



M-kMk)) 



i\k 2 



dp 



(27T) 



z (0(p)0(k-p)0(-k))+AW, (5.5) 



Similarly, multiplying with two fields <fi at different 
wavevectors but equal times, we obtain a relation involv- 
ing the third and fourth moments: 



D (<^(k 1 )0(k 2 )0(k 3 )>^ 



k% 



dp 



i — k\ z 
Po J (2ir)° 

et cycl. 



{ So(k) 
(0(p)0(ki - p)0(k 2 )0(ka)) 



(5.6) 



To obtain this equation we have used the fact that the 
correlation of two fields and one noise r/ vanishes at 
equal time (see appendix). If we proceed in this manner 
we will keep obtaining higher-order correlations in terms 
of the lower ones. We are, however, interested only in 
the lowest-order correlations, so we close the hierarchy 
by replacing the four-point function by a product of two 
two-point functions. Only pairings of p with fc 2 or fc 3 
give nonzero contributions, because of the k\ z prefactor, 
so the dummy index p on the right hand side of ( |5.6| ) 
can only take the values — k 2 and — k3. The two terms 
thus obtained are identical, and we arrive at the following 
relation between the three-point and two-point functions: 

b T.J7k) MkiMk 2 Mks)> = 

1=1 ^ 1 



2zAp 0< 5(k 1 + k 2 + k 3 )So(Wo(fc 2 )So(fc 3 ) S 



kj 



i S (k) 
(5.7) 



Using this in combination with (4.9) we obtain the cor- 
rection to the three-point vertex from sedimentation: 



2i X So(k) 
<5r 3 (ki,k a ,ka) = - s — ^ — 5(k! + k 2 + k 3 ). 



Dpo 2 ^ kr 



t S (k) 



(5.8) 



Combining (5.3) with (5.7), we obtain the correction to 



the structure factor as an integral over wavevectors. In 
the limit of infinite volume, the sum over wavenumbers 
goes from to an upper cutoff determined by the smallest 
length scale in the problem. Our coarse-grained model 
breaks down on scales of order the particle radius R; the 
upper cutoff is therefore Non-dimensionalizing the 
momentum argument in the integral with the — outside, 
i.e. by multiplying all momenta by the average interpar- 
ticle spacing i, we finally obtain 

S(k) = S o (fc)-200 2 S o (fc)^ 

2irt j 

x , R P So (k)S (p)S (\k-p\) 

(27T) 
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x < 



S (p) Sp(fc) S (|p-k|) 
P 2 k 2 (p - k) 2 



S (p) S (k) 5 (|p-k|' 



(5.9) 



where all the wavevectors in the above equation, 
including the fc's multiplying the integral, are non- 
dimensionalizcd by multiplying with I. 

Note that the upper cutoff and hence the integral goes 
to oo in the limit R — ► 0. Let us calculate how the 
structure factor correction behaves in this limit. The 
integrand, for large arguments, ~ | (note that S(p) — ► 1 

for large p), so by power counting, the integral ~ R^^ d K 
This apparent divergence is however controlled by the 
powers of R in the Peclet number, which was R 2d . We are 
thus left with a correction that ~ R( 1+d \ which vanishes 
as R — > for all d. This is consistent with the fact that 
in this limit there are no hydrodynamic interactions. 

The results of this section can also be obtained using 
diagrammatic perturbation theory |pl[ . The correlation 
function calculated to one loop without self consistency 
and the three point vertex at tree level give results iden- 
tical to those obtained above. 



VI. THE STRUCTURE FACTOR 



to be 20, corresponding to 4 ~ 3. The integral has re- 
flection symmetry about 9 = tt/2, so that the structure 
factor at an angle a from the horizontal (6 = 7r/2 — a) 
is the same as that at an angle ~ct from the horizontal 
(9 = 7r/2 + a). This follows from the invariance of the 
integral under k z — > — k z . 

The numerical results are presented in Figs. |l| and ||. 
The first shows the fractional change in the structure fac- 
tor per unit Peclet number as a function of wavenumber 
and angle. The reduction happens to be a maximum at 
the peak itself, and is a monotonic function of the an- 
gle between the horizontal (where it is maximum) and 
the vertical (where it goes to zero). In Fig. |] we show 
the numerically computed structure factor in two dimen- 
sions. The equilibrium structure factor has been taken 
from Monte Carlo simulations of a two-dimensional sys- 
tem of particles interacting via a screened Coulomb po- 
tential, just above freezing p4j . The peak of the equilib- 
rium S(k) was at k = 7.3 (which is close to the average 
interparticle separation fco ~ 2n in the units we are us- 
ing) and the peak value was So(ko) = 4.947. We have 
chosen the Peclet number ^ = 1 for these plots. This is 
a rather large number, and our calculation, which goes 
only upto leading order in perturbation theory is proba- 
bly not accurate for such large driving. The large value 
has been chosen to make the deviations from the equilib- 
rium structure easily visible. 




2 4 6 

The correction 



10 12 

S(fc)-SoO) 
SoCfc) 



16 18 20 

as a function of k for 



different angles with respect to the sedimentation direction, 
in two dimensions: 9 = 90° (dotted line) showing no change; 
9 = 60° (dash-dot); 9 = 30° (dashed line); and 9 = 0° (con- 
tinuous line). Clearly, the correlations are reduced near the 
peak, and the reduction is the largest along the sedimentation 
direction. 



We computed the structure factor of the sedimcnting 
s usp ension by numerically computing the integral in Eq. 
(5.9) fS2[ . We used a 40-point Gaussian quadrature rou- 
tine p3| and took the large k cutoff of the integral 




27m 
R 



2 4 6 8 10 12 14 16 18 20 

FIG. 2. S(k) for a sedimenting colloidal suspension in two 
dimensions as a function of k, at Pe = 1 for four different 
angles with respect to the sedimentation direction: 9 = 90° 
(dotted line) which is the same as the undistorted structure 
factor; 9 = 60° (dash-dot); 9 = 30° (dashed line); and 6 = 0° 
(continuous line). 

It is important to note that the split-peak in S(k) is a 
high Peclet number effect. Since the correction to S(k) 
is peaked at the same place as S(k) itself, there can be 
no shift in the position of the peak of the extremum. At 
low Peclet number, therefore, the structure factor has 
a global maximum at k ~ 7.3, which becomes a local 
minimum at a sufficiently large Pe. At what value of Pe 
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this happens depends on the angle, the required driving 
going to infinity as the horizontal direction is approached. 

We have studie d the freezing transition only in two di- 
mensions (section VII ), but we present structure factor 
data for three dimensions as we ll f Figs. ^ and |^). In 
three dimensions the integral in (5J)), and also the struc- 
ture factor correction, have an azimuthal symmetry. In 
addition there is the reflection symmetry about 6 = n/2 
as in the two-dimensional case. For this case we used a 
30-point Gaussian quadrature routine [p3[ , and checked 
that the answer does not change significantly between 
28-point and 30-point. The structure factor after the 
sedimentation correction is shown in Fig. ^. For the 
equilibrium structure factor we have used Monte Carlo 
results for a \ liquid fl35|| , made continuous by suitable 
interpolation. The Peclet number is once again Pe — 1. 
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FIG. 3. The correction 



10 12 

S(k)-S (k) 
So (*) 



14 16 18 20 

as a function of k for 



different angles with respect to the sedimentation direction, 
in three dimensions; 9 = 90° (dotted line) showing no change; 
9 = 60° (dash-dot); 6 = 30° (dashed line); and 9 = 0° (con- 
tinuous line) . As in the two dimensional case, the correlations 
are reduced near the peak, and the reduction is the largest 
along the sedimentation direction. 
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FIG. 4. The three dimensional S(k) as a function of k 
after sedimentation correction, at Pe = 1 for different angles 
with respect to the sedimentation direction: 9 — 90° (dotted 
line) = the undistorted structure factor; 9 — 60° (dash-dot); 
9 — 30° (dashed line); and 9 — 0° (continuous line). 



VII. FREEZING 

The effect of the corrected three-point vertex as well 
as the structure factor correction have to be considered 
while describing the liquid-solid transition in the pres- 
ence of sedimentation. The latter has a relatively simple 
effect: since S(k) is reduced in all directions, there is 



a tendency towards melting. The correction (5.8) to the 
three-point vertex, being imaginary, has more interesting 
implications. Combining ( |5.8| ) with (4.10b) we conclude 
that the noncquilibrium effects shift the cubic vertex, 



B 



B 



J 1c- 



Po \D 



(7.1) 



E 



So (AO 



where we have scaled all wavevectors by multiplying with 
t he av erage interparticle separation I. With reference to 
( p. 11 ), we can find the cubic term in the Landau expan- 
sion for the nonequilibrium system, 



E 



D 



(7.2) 



E 



i Sa(k) 



Introducing a real dimensionless quantity which is first 
order in Peclet number, 



, 1 = 2( | 



E 



=i so(k) 



(7.3) 



E 



i So(k) 



we rewrite (7.2) as 



b + i(2pn- 1/2 )p 1 



(7.4) 



The correction to the Landau parameter b is entirely 
imaginary, and makes the effective b complex. The free 
energy is a real , positive quantity, so a complex coef- 
ficient B in 2.8 would imply that the complex Fourier 



amplitudes </>(k) pick up phases. In equilibrium, when 
B is real (and positive by construction), the phases of 
</>(— k), (j)(p) and ^(k — p) add up to 27r so as to make the 
cubic term as negative as possible. In other words, if we 
write the amplitude as 







0(k) = |0(k)|e ie ( k) 



then 



0(ki) 
2tt 



0(k 2 )+0(k 3 ) 



(7.5) 



(7.6) 



for sets of vectors that form triangles, i.e. satisfy 
ki + k 2 + k3 = 0. Note that the cubic term is always 
real, whatever the phases of the wavevectors, since for 
any set {ki,k2,ks} there is another with all wavevec- 
tors reversed which has exactly the opposite total phase 
6. However the contribution to the free energy from two 
such terms is 2cosQ, which should be = 2 to minimize the 
free energy. 





FIG. 5. Possible orientations of inner shell reciprocal lat- 
tice vectors in a sedimenting colloidal crystal. The continuous 
lines stand for the distorted crystal, and the dotted lines for 
the original undistorted lattice. The sedimentation direction 
is downwards on the plane of the paper, i.e., antiparallel to 
Gi in the first figure. 



With a complex coefficient as in (7.2), the free energy 
is minimized if the phases satisfy 



0(ki, k2, ka) + arctan 



(2pn- 1 / 2 )(3 1 



2tt. 



(7.7) 



This 'phase shift' in the order parameters may be visu- 
alized in real space as a shift of one set of lattice planes 
relative to its equilibrium position. Note that scattering 
experiments measure only the squared modulus |0(k)| , 
and will not see a phase shift. 

Freezing is determined by the modulus of the modified 
coefficient of the cubic term, which is 



b 2 + 4p 2 n~ 1 /3 1 



(7.1 



The prefactor multiplying (3\ is a number that depends 
on the type of lattice. For a triangular lattice in two di- 
mensions, n — 6 (six reciprocal lattice vectors in the first 
shell) and p = 1 (each reciprocal lattice vector belongs 
to one of the two triangles). Thus 



1,2 + \^ 



(for a planar triangular lattice). 

(7.9) 



Clearly, the magnitude of b can only increase due to sed- 
imentation. Thus, as far as this term is concerned, sed- 
imentation seems to favour ordering. This effect has to 
compete with the reduction of the structure factor and 
only a quantitative analysis can decide their combined 
effect. 

In the rest of this section we present results for the 
simplest case, that of a planar triangular lattice. We 
have described the equilibrium freezing with a one-order 
parameter Landau theory, choosing the Landau coeffi- 
cients to match the structure factor peak height at co- 
existence for our data, i.e. So(k ) = 4.947. Our Landau 
parameters were b — 6, c = 59.364 and, at coexistence, 
a = ciq = 0.20214271, giving a polynomial 



F = 0.10107ar — r 



2.5r 4 



(7.10) 



where r is the reduced order parameter, defined in ( [2.9] ). 

Note that the correction (3\ [eqn. ( [7.3[ )] vanishes if the 
reciprocal lattice vectors remain on the shell |G| = Go, 
and t he st ructure factor remains unchanged. This follows 
from (7.2) and the fact that the sum of ki has to be zero. 
Thus, to benefit from the free energy reduction that an 
increase in b can cause, the lattice may like to distort, 
thus sending one of the vectors out of the shell while 
keeping the sum equal to zero by rotating the other two 
towards or away from each other. We have considered 
two such distortions in our analysis, as shown in Fig. ||, 
which seem to be most likely from the symmetry of the 
problem. It is also imaginable that the lattice distorts 
to an asymmetric state with two degenerate minima of 
the free energy; we have not made a complete scan of 
all possible distortions to rule out such possibilities, but 
looked at only the above two on grounds of plausibility. 

Such distortions apart, the structure factor is itself 
affected by the driving, so (3\ is nonzero even without 
pushing the reciprocal lattice vectors out of the shell. 
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The undistorted crystal is therefore also a strong can- 
didate. The crucial issue is the relative importance of 
the two effects: f3\ and S(k). Though the correction to 
the three-point vertex is first order in Pe, the correction 
to b for small Pe goes as Pe 2 as can be s een simply by 
binomially expanding the square-root in (7.9). For this 
reason, combined with the fact that the numerical value 
of (3% is rather small at modest driving rates, the T3 effect 
is negligible compared to the structure factor reduction. 
For the Landau potential chosen, we have found numer- 
ically that the structure factor reduction makes the free 
energy function go unstable well before b can change sig- 
nificantly. 



0.01 



0.008 



|Gi| = Gq(1 + A), 



(7.11) 
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FIG. 6. The distortion A as a function of Peclet number. 

The physics, therefore, is dominated entirely by the 
structure factor. This implies that sedimentation will 
neccessarily lead to a shift of the liquid-solid phase 
boundary towards the solid side. In other words, since 
the correlations are reduced, a system which was at 
liquid-solid coexistence in the absence of sedimentation 
will choose the liquid phase if sedimentation is switched 
on. 

The reduction of the structure factor peak height is 
maximum along the sedimentation direction and falls 
monotonically as the angle with respect to the horizon- 
tal is decreased, going to zero in the horizontal direction. 
It is therefore expected that one of the wavevectors will 
point along the x direction. This immediately rules out 
the possibility shown in Fig. ||(a) in favour of that shown 
in Fig. ||(b). Even if there is no distortion, this orien- 
tation is the best among the undistorted ones. We now 
ask: is there any reason the lattice may wish to distort 
due to structure factor effects alone? 

Consider the distortion sketched in Fig. ||(b) - the 
reciprocal lattice vector G\ growing slightly larger, 



and the others rotating towards the x axis so as to keep 
the sum of the vectors zero. If the Peclet number is not 
too large, S(k) continues to have a peak at fc even in 
the presence of sedimentation (section VI), so that there 
is a free energy cost to pay for the expanding G\ which 
is moving off the peak. The other two, however, have 
moved closer to the horizontal, i.e. in a direction such 
that the driving affects them less, and the peak height at 
these wavevectors is higher than it would have been had 
they not rotated at all. There is a competetion between 
these two effects, and its not a priori clear which will 
win, so we computed the free energy of distorted crystals 
at coexistence to find the most favourable value of the 
distortion or angle. Note that in any case the liquid is 
the favoured phase, but we can still search among all the 
crystalline phases and find which of them has the lowest 
free energy. This, as we demonstrate shortly, gives us 
some hint as to what to expect at other temperatures. 

We found that for low Peclet numbers the undistorted 
crystal is favoured, but at around Pe — 0.72 a distorted 
crystal becomes favourable (Fig. ||). There seems to be a 
first order jump in what we might call the distortion order 
parameters A. This is not surprising, as we may regard 
the nonequilibrium free energy coming out of our analysis 
as a function of A, whose power series expansion will have 
odd order terms since A — > —A is not a symmetry of 
the problem. 

The distorted crystal appears at rather high Peclet 
number (Fig. |^), certainly well beyond the range of valid- 
ity of our perturbation theory. The reason for the sudden 
distortion is probably related to the original peak ko be- 
coming a minimum and two peaks developing on cither 
side (Fig. |J). The several plateaus that appear at even 
higher Pe are even more intriguing. The intriguing re- 
sults regarding split peaks (Figs. pH) and distortions are 
however somewhat speculative in the absence of a reliable 
theory for large driving. 

We have also compared the configurations shown in 
Figs ||(a) , and ||(b) and found, as expected, that the lat- 
ter has a lower free energy at all values of the driving. 
Thus one immediate consequence of our theory is that 
the crystal, when it forms, will form with the orientation 
shown in Fig. |5|(b), i.e. with lattice planes parallel to 
the sedimentation direction. This agrees with intuition, 
since there are shocks [^4 36,37]] travelling along the sed- 
imentation direction. 
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FIG. 7. Nonequilibrium phase diagram of the sedimenting 
colloidal suspension. The 'a' on the abcissa is the inverse 
structure factor height in the liquid phase of the suspension. 

Finally, one would like to look at the system below 
the equilibrium freezing temperature T c , and find how 
the melting Peclet number depends on T — T c . Strictly 
speaking this requires structure factor d ata in the super- 
cooled liquid, which one should use in (5.9) to compute 
the corrected S(k). Such data, however, is not readily 
available, and in particular for the case of the planar tri- 
angular lattice that we are considering, it is very hard to 
sustain a supercooled liquid for too long in a simulation 
with monodisperse systems. 

We therefore attempt the following approximate the- 
ory. Near T C) the structure factor is that at coexistence 
plus a perturbative correction. The critical Peclet num- 
ber Pe c at which the crystal melts is an object which is 
itself perturbative in T — T c . In calculating Pe c , if we 
neglect the perturbative corrections to So(k), we will be 
making an error only at a higer order. We can therefore 
legitimately do a Landau theory at a T below but close 
to T c using the structure factor at coexistence. The S(k) 
corrections will all be calculated using this So(k; T c ), the 
temperature entering only in the Landau parameter a. 
The result is shown in Fig. ^. 



VIII. CONCLUDING REMARKS 

We have shown that a uniform drift, as in steady-state 
sedimentation, reduces translational correlations in a sus- 
pension of interacting Brownian particles, as a result of 
local hydrodynamic interactions that reduce the mobility 
of a dense region. This reduction is largest in the direc- 
tion of sedimentation. At sufficiently high Peclet num- 
ber, the first peak of the structure factor splits due to 
the fact that the reduction is strongest where the initial 



correlation was largest. 

For the freezing transition we find that the phase 
boundary in a charge stabilized colloidal suspension shifts 
towards the crystal side. The shape of the phase bound- 
ary is quadratic in Pe as shown in Fig. [?]. 

Let us now make some estimates for Peclet numbers in 
realistic experimental situations to see whether the values 
of Pe in Fig. fj] are accessible. The Peclet number for a 
dense suspen sion, at the scale of the particle radius, is 
given by Eq. pT) , which can be rewritten as 



Pe 



Do 
D 



mgR 



(8.1) 



I'o 



provided the hindered settling speed (or the mobility) 
is roughly linear in the volume fraction in the range of 
interest. Taking the specific case of polystyrene spheres 
of diameter 15.5/xto in water Jl5j , we note that this is true 
for volume fraction between 0.1 and 0.3 |58j. The volume 
fraction used in our calculations (-^ *~ 3, from the cutoff 
used for the integrals, i.e. volume fraction — ~ 0.16) 
falls within this range. ^From the data p rese nted in [p8[ , 
the term within the square brackets in (8.1) is found to 
be of order 1.5. 



The first term 



D 



1) is the ratio of the bare 



to the collective diffusion constant, which has a value of 
about 10 near the freezing transition of a colloidal sus- 
pension |39| . The ratio of the interparticle spacing to the 
radius of the polyb alls, (ji), is ~ 3. The product of the 
first three factors in ^j] is thus 10 x 3 x 1.5 = 45. Finally, 
we have to estimate the bare Peclet number ^2^. This 

H — b 

turns out to be ~ 0.5 for polystyrene [|40J. The Peclet 
number for the interacting suspension is therefore 22.5. 
This is rather large, which means that our lowest-order 
theory is probably not reliable here. It does however 
suggest that the effect of steady-state sedimentation on 
freezing should be large for this system. 

In general, in choosing a system on which to test our 
results, two competing needs must be kept in mind. 
On the one hand, the ordering length scale should be 
large compared to the particle size, so that the simple 
concentration-dependence of the mobility (see the re- 
marks towards the end of section [II) in our model is 
valid. On the other hand, the particle size itself should 
be large enough to give an appreciable Peclet number. 

A comparison with shear induced melting is tempt- 
ing, since a perturbation theory similar in spirit to the 
present approach is known to give a shift in the phase 
boundary that is quadratic in the shear rate. There are 
essential differences however: firstly the only significant 
effect of shear is to suppress correlations, whereas for a 
sedimenting suspension, there are nontrivial corrections 
to the three point vertex. The effect of these latter correc- 
tions happens to be negligible on the freezing transition 
as we have shown, making the phase boundary shift in a 
manner similar to the shear case. Secondly, the sheared 
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colloidal crystal at low shear-rates is a slightly distorted 
crystal essentially stuck in a statically sheared configura- 
tion. At low shear therefore (where perturbation theory 
is valid) , the crystal is essentially in an equilibrium state, 
and its properties can be derived from an cqulibrium pic- 
ture |J. For uniform driving, however, there are nonzero 
currents even at low Peclet numbers, making it a more 
genuinely nonequilibrium problem. 

An interesting feature of our theory is that the sign of 
^ does not figure in the results. Thus even if we have a 
suspension whose settling speed increases due to many 
body effects, our theory will still predict a reduction 
in the correlations and a consequent shift in the phase 
boundary towards the crystal side. 

Finally, a few caveats. In ||l|] we have shown that a 
scdimenting colloidal crystalline lattice is likely to be un- 
stable at long wavelength. The theory proposed in this 
paper is probably nonetheless applicable to the formation 
of the crystalline phase at least at low Peclet numbers, 
since the physics of freezing is centred primarily around 
the ordering wavenumber. The coupling of concentration 
and shear (tilt) considered in 41 4^] may also be impor- 
tant in the liquid state. The detailed consequences will 
be different because shear stresses have a finite lifetime 
in the liquid, in contrast to that in a solid where it lasts 
for ever if the strain is present. The effect of the sedimen- 
tation instability is thus very different in a liquid. This 
is an additional aspect that we have totally ignored in 
the above treatment, and one that a more sophisticated 
approach must address. However, further development 
of the theory should be motivated by experiments, which 
we eagerly await. 



APPENDIX A: FIELD-NOISE CORRELATIONS 
AT EQUAL TIME 

Consider a set of Langevin equations of the form 



<t>a = f(<t>) +ria, 

(Vc(t)V0(t')) = 2k B TT aP 5(t - t'). 



(Ala) 
(Alb) 



where <fi stands for {(pi, 4>2, ■ ■ • and / can be any 
arbitrary function, in general nonlinear. We shall calcu- 
late the field-noise correlations using Novikov's theorem 
[fl3| flsf , which states that for any functional of the noise 
AM, 



(Av a (t)) = 2TT a0 ^ 



(A2) 



We first calculate the average of the product of one field 
with one noise, which by (A2) is 



(A3) 



0a(*)= f dt'f (#'))+ / dt' Vx (t') 

Jt J to 



(A4) 



where the rj\ dependence in the first term is present 
through 4>(t), which depends on the noise amplitudes at 
all past times but not those of the future. The derivative 



Sr) X (V) 



(A5) 



vanishes for t < t' by causality. To calculate the equal 
time value of this derivative, let us calculate it first for 
t = t' + e for infinitesimal positive e. For t > t' , 



SvxW) J to 9 h <%(f) 



(A6) 



For t — t' , all the arguments t" in the integral are times 
before t', so the functional derivative inside the integral, 
and hence the integral itself vanishes. Thus the forward 
derivative is simply equal to 1. The backward derivative 
is zero, therefore 



6<p a {t) 
5r, x V 



6{t-t'). 



(A7) 



At equal time this is discontinuous, but we recall that 
the delta function noises are nothing but limits of Gaus- 
sians, so the theta function is not really a discontinuous 
function at the origin, but rather rises smoothly from 
to 1 taking the midpoint value ^ at the origin. Thus 



<%(*) 



which, combined with (A2) gives 



(4> a (t)r)\(t')} = k B TY aX . 



(A8) 



(A9) 



This can be shown with more rigour by starting from a 
Gaussian distribution, calculating the equal time corre- 
lation and then taking the width of the Gaussian to zero 
|Q . In section [v| we have used an analogous result for 
reciprocal space variables. 

Let us now consi der the product of two fields and a 
noise. Again, using ( |A2| ), we have 



Svx(t) 



(A10) 



The equation ( |Ala|) can be formally solved to give 



where we have used the chain rule for functional deriva- 
tives. Now each of the partial derivatives = ^, therefore 

(Mt)Mt)vx(t)) = \ (Mt) + Mt)) ■ (ah) 

In section [v] we do a perturbation theory about the liquid 
phase, calculating liquid state correlations, so the first 
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moments of <fi are an zero, thus making all averages of 
the form (4>4>rj) vanish. 

For completeness we explicitly check the vanishing of 
t he first moment in the presence of driving. Averaging 
(4.4) over noise, and setting the time derivative of the 
field average to zero (steady state), we have 



Dk 2 
S (k) 



i—k z 



dp 

(2^ 



<0(p)0(k - p)} . (A12) 



The second moment on the right is nonzero only if — p = 
k—p, i.e. k = 0, so that the first moment (</>(k)) vanishes 
for all k ^ 0. But since 4> is a deviation variable, its k = 
component is zero by construction. 
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